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We prove existence of a shape and boundary condition independent thermodynamic limit for 
fluids and solids of identical particles with electric or magnetic dipole moments. Our result applies 
' to fluids of hard core particles, to dipolar soft spheres and Stockmayer fluids, to disordered solid 

(£) , composites, and to regular crystal lattices. In addition to their permanent dipole moments, particles 

^"*^ ■ may further polarize each other. Classical and quantum models are treated. Shape independence 

depends on the reduction in free energy accomplished by domain formation, so our proof applies only 
in the case of zero applied field. Existence of a thermodynamic limit implies texture formation in 
spontaneously magnetized liquids and disordered solids analogous to domain formation in crystalline 
solids. 
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"■^J , Thermodynamics normally assumes a free energy density F/V exists and is independent of system volume V 
£^ ■ and shape. Verification of these properties is impeded by the explicit dependence of the partition function Z = 
exp (—F/JcbT) on these very quantities. Ruelle 0] and Fisher p| proved the existence of thermodynamic limits for a 
i large class of fluids and solids with interactions that fall off faster than r~ 3 at large separation. For such systems the 

free energy contains a boundary independent, extensive (proportional to system volume) component and a boundary 
dependent, sub-extensive (less than proportional to system volume) remainder. Consequently, in the limit of infinite 
volume the free energy density approaches a finite, boundary independent, limit. 

The interaction energy between dipoles falls off precisely as r -3 , seriously complicating the thermodynamic limit. 
Volume integrals of this interaction (required to calculate the total interaction energy H) converge only conditionally 
because the power of r with which the interaction decays matches the dimensionality of space. This paper considers 
systems with electric or magnetic dipole interactions. The electric and magnetic cases resemble each other closely. 
For convenience we carry out our discussion in the context of magnetism, then address electric analogues near the 
. end of the paper. 

Long-ranged dipole interactions may create shape dependent internal demagnetizing fields that increase the system 
free energy. Boundary conditions on the surfaces may influence the strength of these demagnetizing fields. The 
reduction in demagnetization energy when uniformly magnetized regions break into smaller domains is the key to the 
very existence of a thermodynamic limit in zero magnetic field. Griffiths || used the reversal of magnetization in a 
domain to prove existence of a thermodynamic limit independent of shape for dipolar lattices. We generalize that 
proof to include fluids and disordered solids. Certain conditions are required on the "residual interaction" H R , defined 
as the total interaction energy TL minus the magnetic interaction energy 7i M . Our proof, like Griffiths' original one, 
is valid only for zero applied field because of its reliance on magnetization reversal in domains. 

The following section of our paper describes the or igin of demagnetizing fields, leading to a shape dependent free 



energy in the presence of an applied field. In section II A we conjecture a simple functional form of the free energy, 



X. 

absorbing all shape dependence into a de mag netizing energy, implying a conventional thermodynamic limit for the 



remaining part of the free energy. Section II B describes how the system achieves a thermodynamic limit in zero field. 



Then, in section III, we outline the formal thermodynamic limit proof, which relies o n upp er and lo wer b ounds on 



the free energy. We illustrate these bounds for stable and tempered systems in section III B. Section [II C discusses 
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the difficulty dipolar interactions cause due to their lack of tempering, and how that difficulty may be overcome. 
Section IV extends the proof to a variety of interesting specific models, starting with identical hard core particles, 
then treating dipolar soft spheres, the Stockmayer fluid and polarizable particles. We treat both classical and quantum 
versions of all these models. Section ^ addresses the analogous problems for electric dipoles. Finally, in section Jv| we 
summarize our results, discuss some observations about the implications of a thermodynamic limit for spontaneously 
polarized liquids, and conclude with some interesting dipolar systems which lack a thermodynamic limit. 



II. DEMAGNETIZING FIELDS 



Long range dipole interactions create demagnetizing fields that cause the shape and boundary condition dependence 
of free energy in the presence of an applied field. In zero applied field the demagnetizing field must be handled with 
care to prove the shape and boundary condition independence of free energy. This section describes qualitatively why 
materials lack a shape independent thermodynamic limit in an applied field, and how the limit is restored in the 
absence of applied fields. We conjecture a modified form of thermodynamic limit that may hold in an applied field, 
with all the shape dependence restricted to an effective internal field plus an explicit demagnetizing energy term. 

Volume and surface distributions of magnetic poles cause demagnetizing fields. Consider a sample of magnetic 
material contained in a region of space of volume V. Let M(r) be the spatially varying magnetization in the sample. 
Poles arise at surfaces wherever the magnetization has a component normal to the surface S, as indicated in figures |a 
and [j]b. A magnetization with non-zero divergence produces a charge density in the bulk. The demagnetizing field 
takes the form 

U D (r) = [ dV (M(r') • n(r'))-^— ^ - / dV (V • M(r')) J— - ^ (1) 



where fi(r') is the outward normal at any point on the surface. The spatial arrangement of surface charges depends 
on sample shape, and the divergence of M(r) in the interior depends on the magnetization texture, so H^i(r) is a 
function of sample shape and magnetization texture. The demagnetizing energy 



~ J d 3 r H D (r) • M(r) = i- f d 3 r \H D (r)\ 2 > 



E D = --\ d 3 r H D (r) • M(r) = — / d 3 r \U D (r)\ 2 > (2) 

depends explicitly on the shape and magnetization texture of the system through the demagnetizing field H.o(r). 

In the special case of magnetization uniform throughout the sample, the demagnetizing field H£>(r) comes only 
from the surface, because the divergence term in equation (|) vanishes. However, Ho(r) does not vanish as volume 
increases at fixed shape. This is because the I jr 2 fall-off of the field from each surface charge is exactly offset by the 
r 2 growth of surface area, and hence the number of surface charges. As a result, Hd(i*) is independent of the volume 
and the demagnetizing energy Ed is extensive. 

For the special case of a uniformly magnetized ellipsoid, Ho is constant within the ellipsoid and equals 

H D = -4ttD • M, (3) 

where the tensor D is the demagnetizing factor of the ellipsoid [Q. D is non-negative definite, and its trace equals 
1 . When the magnetization lies along a principal axis of the ellipsoid, D is simply replaced by one of its eigenvalues 
< D < 1. For a magnetization parallel to a highly elongated needle shape, the demagnetizing factor D — because 
the surface poles appear only on the tips which are small and far removed from the bulk. Another special limit is that 
of magnetization normal to a flat pancake shape. This yields the maximum demagnetizing effect, since the surface 
poles appear on a large surface close to the bulk, so D = 1 in this case. 



A. Shape dependence in a field 

When a system is placed in an external field Ho, surface poles arise because the internal magnetization tends to 
align with the applied field. There are two important contributions to the resulting shape dependence of the free 
energy. One is the explicit shape dependent energy (||), the other is due to the shape dependence of the internal field 

H= H + H D . (4) 

For highly elongated sample shapes, in the absence of demagnetizing effects, we expect a thermodynamic limit for 
the free energy. Define /j ni (Ho) as the free energy per unit volume of a system in the limit as total volume V —>■ oo. 
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The limit must be taken within ellipsoidal shapes for which the length parallel to the field Ho grows faster than the 
orthogonal directions. Because there are no demagnetizing effects present, we call this free energy density the intrinsic 
free energy density in a field. 

For more general shapes Hd is non-zero, and may vary in space. We conjecture that the shape dependent free 
energy F s h ape may be expressed in terms of the intrinsic free energy density as 

Eshape (Ho) = J d 3 r /i„t(H(r)) + E D (5) 

up to corrections that grow less rapidly than the volume. Equation (Q) gives the internal field H and equation (|^) 
gives E D . 

Shape dependence of the free energy implies shape dependence of the measured paramagnetic susceptibility. Assume 
that the magnetization M is related to the internal field H by an intrinsic (volume and shape independent) linear 
susceptibility Xint according to 

M = Xint H. (6) 

Consider applying an external field Ho parallel to a principal axis of an ellipsoidal sample. Because of the demagne- 
tizing effect, the internal field H is weaker than the applied field. Eliminating and H between equations (^|), (Q), 
and (^) yields the shape dependent measured susceptibility 

M = Xshape H (7) 

where 

Xint /„% 

Xshape — Z j Z — ■ \o) 

Note that the measured Xshape has a maximum value equal to Xint when Ho is parallel to the long axis of a highly 
prolate needle-shaped ellipsoid. For any other geometry the demagnetizing effect reduces the measured susceptibility. 



B. Shape independence in zero field 

Now consider a ferromagnetic material in zero applied field. If the magnetization were constant (Fig. |l]a), surface 
poles would create shape dependent demagnetizing fields and raise the energy as described in equation (||) . A uniformly 
magnetized body lacks a shape independent thermodynamic limit! 

Alternative magnetization configurations reduce the demagnetizing energy. One possibility (Fig. |f|b ^reverses mag- 
netization in subregions so that the fields from surface poles tend to cancel. Another possibility (Fig. [lp) rotates the 
magnetization so that it is always tangent to the surface. In each case, the reduction in energy is proportional to the 
system volume L 3 , where L is a typical linear dimension. The energy increase arising at a sharp domain wall (Fig. |l]b) 
should be proportional to the domain wall area L 2 . The energy of a vortex line (Fig. |f|c) should be proportional to 
Llog(L/a), with a related to the vortex core size. The magnetization texture may avoid a vortex by escaping into 
the the third dimension B near the core. Such textures contain either surface poles f| or point defects Q. For 
sufficiently large systems, the extensive L 3 reduction in the demagnetizing energy dominates the sub-extensive defect 
energies, and domain wall, vortex, surface pole or point defect formation is favored in that it lowers the free energy. 

As we take the V going to infinity limit, the demagnetizing energy density Ed /V approaches zero for the most 
favorable magnetization, i.e., the one which minimizes the energy. Such a magnetization permits a shape indepen- 
dent thermodynamic limit for a ferromagnet. The free energy density for an arbitrary shape with its nonuniform 
equilibrium magnetization texture equals the free energy density of a highly elongated needle-shaped ellipsoid with 
uniform magnetization parallel to the long axis, because D tends to zero for the needle shape. When calculating 
magnetic energies or free energies it may be convenient to impose the needle-shape and assume uniform polarization. 
Alternatively, "tin foil" boundary conditions || may be used to neutralize the surface poles. 

Why is zero applied magnetic field essential for a thermodynamic limit? In an applied magnetic field the energy 
cost for flipping a domain becomes proportional to the domain volume (and grows proportionally to L 3 ) rather 
than the smaller domain wall, vortex or other defect energy. The most favorable magnetization texture now has a 
demagnetizing field, and the free energy re-acquires a shape dependence. 

The above discussion explains how domain formation removes the demagnetizing energy density Ed /V, permitting a 
shape independent free energy for zero field fcrromagnets. This argument does not apply to the zero field susceptibility 
of a paramagnet. The shape dependent susceptibility Xshape governs fluctuations in the average magnetization of the 
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entire sample. When this average fluctuates from zero, demagnetizing fields increase the free energy and oppose the 
fluctuation. Reduced fluctuations imply a reduced susceptibility that depends explicitly on shape through equation (||). 

Still, we expect shape independent values of magnetic permeability p, = 1 + ^irxint (or dielectric constant e in the 
case of electric polarizability). This can be understood by expressing the permeability in terms of spatial integrals of 
correlation functions p|JTo|]. These correlation functions contain short-ranged, shape independent, components, and 
long-ranged, shape dependent components. The permeability depends only on the short-ranged, shape independent, 
part of the correlation functions. 



III. PROOF OF THE THERMODYNAMIC LIMIT 



This section explains how we prove thermodynamic limits. First, we state required bounds on the free energy and 
explain how these bounds are used to prove the existence of a thermodynamic limit. Then, we show how to prove the 
necessary bounds on the free energy for classical systems which are stable and tempered. The se sections are rather 
brief and formal, and simply review methods introduced previously . Then, in section III C we show how to treat 
systems which include unstable and non-tempered dipole interactions. 



A. Conditions on the Free Energy 



Consider an N particle system contained in a region 1Z of volume V. Taking the thermodynamic limit for the free 
energy means constructing a sequence of sufficiently regular regions jl , with increasing volume, so that the number 
of particles N divided by V approaches a definite value p as the volume V tends to infinity. A limit is said to exist 
for the free energy density if the free energy F divided by the system volume V approaches a limiting value / as the 
volume tends to infinity. The requirement of regularity B prevents the regions 1Z from getting too thin or constricted. 
We also introduce a model-dependent density p c that ensures the particles can fit into the available volume when 
N/V is less than p c for sufficiently large finite N. 

Two conditions on F suffice to prove the thermodynamic limit. 

1) The free energy F should satisfy the lower bound 

F>Vf L (N/V) (9) 

for N/V < p c where Jl(p) is some finite valued function. 

2) Consider a system composed of two subsystems, 1 and 2, containing N% and N 2 particles, respectively. The 
particles in subsystem 1 are confined in a region IZi with volume V± and those in subsystem 2 are confined in a region 
1Z 2 with volume V^. The two regions 72. i and 7^2 are separated by a distance of at least d from each other. Provided 
that d > c?q , for some fixed distance do, the free energy of the system should satisfy an upper bound 

F<F 1 +F 2 + A 12 , (10) 

where F\ and F 2 are the free energies of subsystems 1 and 2 in isolation and Jll[] 

A 12 = (N^N^tUB/d 3 ^, (11) 

with constants lob < oo and e > 0. 

These bounds suffice for proving the existence and shape independence of the thermodynamic limit. Break an 



arbitrarily shaped system into many smaller subsystems. The upper bound (10) bounds the total free energy in terms 
of the subsystem free energies. Because the upper bound applies regardless of the relative positions of the subsystems, 
provided d > d Q , the original system shape does not enter this bound on total free energy. The lower bound (||) 
guarantees that the free energy density F/V reaches a finite limit as the total volume V — > oo with N/V — > p. 
Because Fisher Q explains this method in great generality, we need not reproduce his effort here. 

B. Classical Stable and Tempered Systems 

The free energy of a classical system of N identical particles in a volume V is F = — ksT log Z, where Z is the 
partition function 
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The energy TL is a function of particle center of mass positions , and internal coordinates tti . In the expression for 
Z, is the integral of d£li over all its possible values. Internal coordinates depend on the type of particle and may 
include orientation of the particle and direction of magnetization. For solids, particle center of mass positions and 
particle orientations are fixed, and the principal remaining variable is direction of magnetization. 

We distinguish between two types of particle: superparamagnetic particles Jl2[ , for which the direction of magne- 
tization rotates independently of the particle axes; normal particles, for which the direction of magnetization is fixed 
relative to the particle axes. For normal particles, we do not include direction of magnetization as an independent 
internal variable, because it is a function of particle orientation. In practice, superparamagnetic particles exhibit 
a "blocking temperature" below which the direction of magnetization becomes locked to the particle axes, and the 
particles become normal. In the specific models discussed below, we assume we are below the blocking temperature 
except where we explicitly invoke superparamagnetism. 

Note the explicit dependence of Z on the system shape through the limits of integration f or the in equation (Mj 



The free energy F inherits this shape dependence. The conditions on F stated in section [III A| guarantee that the 
shape dependence is contained entirely in a sub-extensive term. Achieving the desired lower and upper bounds on 
free energy depends on properties of the interaction energy TL. This section describes sufficient conditions to prove 
each bound. 

The lower bound @ holds for potentials that are stable in the sense that 

TL > -Nlu a (13) 

with uja < oo a constant. Just substitute the lower bound ( |l3| ) for TL into the partition function ( |l2| ) to obtain the 
lower bound (||) on F, with the function 

f L (p) = pk B T log p - puj A . (14) 

To prove the upper bound ([To| ) , consider the interaction of two subsystems separated by distance d as described in 
section III A . Write the total energy TL in the form 

H = H!+H2 + H 12 , (15) 

where TLi and TL2 denote the energies of each system by itself, and TL12 is the interaction energy between the two 
subsystems. The upper bound holds if the interaction TL12 satisfies the weak tempering condition [|IlJ 

W12 < A 12 (16) 

with A12 as defined in equation (|Tl|), for d larger than some constant do. Substitute A12 for TL12 in the total interaction 
energy ([l5|) and evaluate the partition function ( |l2| ) to derive the upper bound ( |io| ) on F. 



C. Dipolar Systems 

The remainder of this paper considers systems whose Hamiltonians include dipolar interactions in addition to stable, 
tempered interactions of the type described above. The dipole interaction, by itself, is neither stable nor tempered. In 
this section we explain how additional repulsive interactions may stabilize the system, and how the upper bound ( |l0| ) 
may be proven d espi te the lack of tempering. The ideas introduced here are applied to a wide variety of specific 



models in section IV 



Split the interaction energy TL into two components: 

TL = TL M + TL R . (17) 

The non-magnetic part of the interaction, TL R , we call the residual interaction. The magnetic interaction between the 
N particles takes the form jj] 

TL M = f / *r [ ^( M W- M f) _ 3(M(r) ■ (r - rQ) (M(r ; ) ■ (r - r')) j 

£sUK 4 I l r - r 'l 3 l r - r 'l 5 J 
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where Vi and Vj are the regions of space occupied by the magnetic material of particle i and j. The magnetization 
distribution of the i particle is M(r) for r inside Vi, and for the j th particle is M(r') for r' inside Vj. Implicitly, H M 
depends on the particle center of mass positions and the particle orientations through Uj and Vj , the regions of space 
occupied by the particles. For superparamagnetic particles, the direction of magnetization is an internal coordinate 
for each particle, while for normal particles, the magnetization is determined by the particle orientation. Thus, H M 
is a function of particle positions {r^} and internal coordinates {f^}. 



For the moment we consider only permanent magnetization (polarizable particles are discussed in section [V C ) , 
and we assume the magnetized volumes of the particles are non-overlapping. The 1/r 3 dependence of the magnetic 
interaction H M violates tempering because of its slow decay at long range, and risks violating stability because of its 
divergence at short range. 

We demand stability of the total interaction H to enforce the lower bound (^) on the free energy. Because of the 
diverging short-range magnetic attraction, we need residual interactions that are sufficiently repulsive at short range 
to overcome the magnetic attraction. Hard-core particles, and soft-core particles with energies that diverge faster 
than 1/r 3 , satisfy this requirement, as we prove later in section IV. 



11111 -L / I . UU' \J 1U1J U111U 1 VV^Lllll^lll^lllI, UIU «V Jr*^- *-* * V> l«-"UVjl 111 UVjVjUHJII 1 Y . 

To achieve the upper bound ( |10|) on the free energy we demand that the residual interaction H R be tempered 
and we exploit symmetries (if present) to handle the non-tempered magnetic interaction H M . Our strategy limits 
our proof to models possessing the required symmetries and tempering of residual interactions. Models lacking these 
characteristics may still possess a thermodynamic limit even though we cannot prove it. The symmetries we require 
are broken by applied magnetic fields. 

Consider two subsystems such that the iVi particles in region IZi are separated by at least a distance d > do > 
from the N 2 particles in region 7£ 2 . Let Hi and H2 be their respective Hamiltonians. Define the interaction energy 
between the two subsystems by 

H12 = H - Hi - H 2 = H12 + H-12- (19) 

Let F(X) be the free energy of the combined system when the Hamiltonian is Hi +H2 plus a scaled interaction XHi2- 
Because F(X) is a concave function (that is, F"(X) < 0), it is bounded above by 

F(X) <F(0)+AF'(0), (20) 

where the right side is a line tangent to the graph of F(X) at A = 0; here F'(X) and F"(X) are the first and second 
derivatives. As a consequence, the free energy F(l) of the fully interacting system satisfies the Gibbs inequality [jl3| 

F(1)<F(0) + (H™+H R ) X = , (21) 

where F(0) = Fi + F 2 is the free energy of the non-interacting subsystems, and the classical ensemble average of any 
quantity Q takes the form 

JVi N 2 

<«« = n^-N^ZA /,, II / h II 8e-<— V-*. (22) 

Because the residual interaction H R is tempered, therefore H R 2 < A12. Its ensemble average, likewise, is bounded 
above: 



(«&>a=o < A 12 . (23) 



We will show that 



( H& >a=o = 0. (24) 

Combining the bound ( pl[ ) with the ensemble averages ( p3] ) and (24) proves the upper bound (|l^) on the fully 
interacting free energy F~(l). 

To establish ( pi|) we employ what we call a 9 operator, a map of the coordinates of a system onto themselves in a one- 
to-one manner satisfying the following conditions. It leaves the center of mass position r,; of each particle unchanged, it 
maps the internal coordinates Qi onto themselves in a way which leaves the integration measure Yii=i unchanged, 
and it leaves the Hamiltonian H invariant. In addition, when a system consists of two subsystems and is applied 
to one but not the other, it reverses the sign of the magnetic interaction Hf% between them. Specific models may or 
may not possess such an operator. When a system is stable and possesses a 6 operator, we can obtain the desired 
upper and lower bounds on F to prove a thermodynamic limit. 

If a 9 operator exists, it can be used to establish ( p4] ) in the following way. Set Q = Hf^ in (|22|), and consider the 
change of variables produced by applying the 9 operator to subsystem 1 but not subsystem 2. This change preserves 
the integration measure, but reverses the sign of the integrand, since Hi and H2 are unaltered, but W^f changes sign. 
Consequently, the integral is equal to its negative, so it is zero. 



G 



IV. MODELS 



Section [II introduced a general strategy for proving thermodynamic limits of permanently magnetized classical 
particles. The following section applies that strategy to a variety of models. We start with identical hard core particles, 
then treat dipolar soft spheres such as Stockmayer fluids. We then modify the proof to cover polarizable particles, 
and then treat quantum systems. Depending on the particular system, the greater challenge may lie in demonstrating 
the lower, or the upper, bound on F. 



A. Identical hard core particles 

Consider a collection of N identical, uniformly magnetized, hard core normal particles of volume v, and fully 
contained within a region of space 1Z of volume V. The magnetization M(r) is constant in magnitude for r in volume 
Vi of particle i and vanishes when r is not inside a particle. Inside particle i the direction of M(r) = M; depends on 
the orientation fli of the particle. We require that the region 1Z have a regular shape (2| and be large enough so that 
all particles fit inside the region without overlapping. Thus, we restrict the number of particles so that the packing 
fraction (f> = Nv/V is not too large. In particular, we assume a packing fraction <fi* > exists for which particles may 
be packed with any < 4> < 4>* into any sufficiently large and regular volume. 

Write the Hamiltonian as 

n = n M + H HC . (25) 

The hard core interaction, 

ruHC _ \ if no particles overlap, 
1 +oo otherwise, 

prevents any overlap between particles. For non-overlapping configurations the magnetic interaction Tt M is as in JT§|). 
For the special case of hard core spheres the expression (|l^) reduces to the simpler form 

i<j=i r « 

where fi i and Yi are the dipole moment and position of the i th particle; /x i is the integral of M(r) over the volume of 
the i th particle; — Yi — Yj, and is the unit vector along r.^ . 

The hard core interaction ( |26|) by itself provides an example of a stable and tempered interaction. Since H HC > 0, 
it obeys the stability condition jl^ ) with uja = 0. Griffiths |3| proved the lower bound 

n M > -^r, (28) 

for non-overlapping dipolar spheres with radius R, and dipole moment [i, regardless of their positions and orientations. 
Because of the hard core repulsion ( p6|) we achieve s tabili ty ( |l3| ) with uja — M 2 /2i? 3 for dipolar hard spheres j2^). 
The lower bound on F follows as discussed in section III B . We now generalize the proof of stability (|l^) , and thus a 



lower bound on F, to particles of all shapes. 

To prove stability we make use of the positivity of field energy. Adding the magnetic self energy of each particle to 
7i M gives the total energy of the whole system, considered as one magnetization distribution, 

N 

H T = H M + E f S = -5 / d 3 Y Hz,(r) • M(r). (29) 
i=i 2 ■> 

Here Hc(r) is the field, due to all particles, defined in equation ([]]) and 

Eff = -~Mi- J rf 3 rH lfl (r), (30) 

where Hj£>(r)is the field from magnetization Mj of particle i with volume v$, obtained by substituting M, for M(r) 
in equation (|l|) and integrating over the surface and volume of particle i. 
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We use equation (^9|) to place a lower bound on 7i M by a method similar to that of Griffiths 
magnetization distribution M(r) and the field Ho(r) caused by it 

~l [d 3 v H d (t) • M(r) = i- fd 3 r |H D (r)| 2 > 0. 



Hence 



N 



For any 
(31) 

(32) 



Brown Q] rewrites the self energy in (^) as 

Ef f = 2-K^DfMfM\ Vl 



(33) 



where Di is the demagnetizing tensor of an "equivalent ellipsoid" ; it exists for a particle of any shape, and k and I 
index the components of and M. Since is positive definite, with trace equal to 1, 



E* elf < 2irM 2 Vl . 

Since all particles are identical, the magnetic interaction satisfies the lower bound 

H M > -2nM 2 Nv. 



(34) 



(35) 

Thus we confirm stability ([l3|). The lower bound (^|) on the free energy follows with to a — 2nM 2 v. 

For proving an upper bound on the free energy, notice that the hard core interaction ( p6| ) is tempered, equation (|l^) , 
with any do > 0. We ide ntify the hard core interaction ( |2lf ) as a residual interaction TV^. The key to our proof of an 
upper bound in Section 1IIC was reversing the sign of the magnetic interaction energy TCf^, without changing Wi, 
by applying an operator 9 on subsystem 1. For particles with permanent magnetization fixed relative to the particle, 
a rotation of each particle can reverse the direction of magnetization. Such a rotation keeps the residual interactions 
unchanged only if the particle shape has an axis of 2-fold symmetry perpendicular to its magnetization. Hence at least 
one operator 9 exists, and our proof applies, for systems of identical particles with the required rotational symmetry 
in shape. 

Some kinds of small particles, including many used in ferrofluids p 4 fl , exhibit superparamagnetism. Dipole moments 
rotate by Neel relaxation [pTs|] , or possibly quantum tunneling ||lq| , without requiring rotation of the particle itself. 
To describe the superparamagnetic classical particles, one includes in fij, in addition to the Euler angles, a discrete 
variable Oi — ±1 specifying that the magnetization is parallel (+1) or opposite (—1) to a direction fixed in the particle, 
and f dQi includes a sum over Oi. The 9 operator is the map Oi — > —Oi, applied to every particle. For a quantum 
particle, the corresponding operation is time revers al applied to the particle's magnetization. In either case, the 9 
operator has the properties specified in section III C , so the argument given there shows that any system of hard core 
superparamagnetic particles, with any shape of particle, has a thermodynamic limit. 



B. Dipolar systems with central forces 



Consider a system of particles interacting with Hamiltonian 

H = H U + H {n) + H central (36) 



where H is the point dipole interaction (07 



N A 



r 

i<j=i 



4-1 'v 



is a repulsive interaction with A > and n > 3, and T-t central is any stable ( |l3| ) and tempered (|T^) potential that 
is spherically symmet ric. D efine Ti^ + Ji central to be the residual interaction TL R . The upper bound ( |l0| ) follows 



exactly as in section IIIC because TiS n ^ and J{ central are both tempered and rotationally invariant. The proof of 
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a lower bound (^J) for such systems is more complicated than for hard core particles because there is no minimum 
distance of separation between the point dipoles. 

To prove stability ( fl3| ) and hence a lower bound (0) , it suffices to prove that H M + ft (n) is stable, since J{ central is 
stable by assumption. Consider some configuration of a finite system of N particles. Let 2Ri be the distance from the 
i th particle to its nearest neighbor. The magnetic interaction energy remains unchanged if each particle is replaced 
by a sphere of radius Ri with uniform magnetization and the same dipole moment \x. The self energy of such a 
sphere is j2R\. Since H M + (/i 2 /2) YltLi ^-1^ > by positivity of field energy (see discussion in section IV A ), and 
H (n) > {A/2)J2? =1 (2Ri)~ n , we write 

N 2 N A 

H M + n (n) > _y_t_ + (38) 



Define a generalized mean 

1 N 1 

G W = (^ElfT)^ (39) 
Using the property that G(t) increases monotonically for positive t JlTj we find 

i=l 1 i=l 1 



for n > 3. Combining equation ( poj ) with equation (^) we write 

H M + H (n) > -7V(^X - ^X"/ 3 ), (41) 

where X = (1/JV) J2iLi R *~ 3 ■ Thc bound @ for U M +H [n) has a minimum because n > 3 and A > 0. In particular 

w ^ + w (n)>_ JV[ ^ ( !^ )( 32V )75 i ST]< (42) 

Z TZ TlJ\ 

Our model ( |36| ) is therefore stable. 

Let's apply this general proof to some special cases. Our proof applies to generalized Lennard- Jones particles with 
dipole interactions. The Hamiltonian for such a system is 

H = H M + H LJ , (43) 

where the Lennard- Jones potential is 

H " = Y — -Y-°-. (44) 

■ , - 1 ii . , . 1 ii 

i<j J ?<j J 

Ruelle showed that generalized Lennard- Jones potentials with B, C > and n > m > 3 are stable. To demonstrate 
stability including thc dipole interaction, divide the repulsive term into two positive pieces, B = A + B' . Attribute 
B' to a new (but still stable) Lennard-Jones potential and use the remainder A to define the repulsive potential 1iS n \ 

H = H M +n in) +H LJ ' . (45) 



The Hamiltonian in fl4q ) is special case of our model (p6[). The Stockmayer fluid |19| is the case with n = 12, m — 6, 
and hence will also have a shape independent thermodynamic limit. Dipolar soft spheres |lg| ] are the trivial case C=0 
and U central =0. 
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C. Polarizable Particles 



Consider a system of identical hard core particles that contain permanent magnetic moments but are further linearly 
(i.e. proportionally to the local field) polarizable. The simplest model for such systems is the dipole-induced-dipole 
(DID) model JlC| , pO| ]. The model consists of spherical particles with a point dipole moment at the center, and with 
the induced polarization an additional point dipole moment of strength aH. This model lacks stability in general. 
For example, with two spherical particles of polarizability a and hard core radius i?, stability is lost when 

The model fails because of its assumed induced point dipole. The point dipole applies rigorously only to an infinitesimal 
volume. However, the polarizability a necessarily vanishes in the limit of zero volume, due to the self-induced 
demagnetizing field. Finite size particles can have a ^ 0. However, the DID model omits multipole moments due 
to the spatial variation of fields and magnetization inside the particles. Higher order multipole interactions between 
particles become important when the particles approach each other J22|,^3| , and are required for stability. 

Consequently, we work with a more physically realistic model [^of . By incorporating the full magnetic interaction 
(dipole and higher moments) and the spatial variation of fields within a particle, our model satisfies stability in 
general. This model represents the polarizability of atoms and molecules more accurately than the DID model. Each 
particle has a permanent magnetization density which is constant in the interior of the particle, but whose direction 
is determined by the orientation of the particle. (For example, imagine that the particles are prolate ellipsoids with 
magnetization along the long axis.) The particles are "hard", so that their volumes cannot overlap. Consequently, 
the permanent magnetization is a vector field M p (r), equal to zero unless r is inside some particle, where it takes on 
a value whose magnitude is independent of the particle but whose direction is tied to the particle's orientation fi. 

In addition, each particle contains linearly polarizable material giving rise to an induced magnetization 

M*(r) = X (r) • H(r), (47) 

where H(r) is the total magnetic field at r, and the susceptibility tensor x( r ) is zero unless r is inside some particle, 
where its value is independent of r but tied to the orientation of the particle: that is the principal values of x an d 
the relationship of the principal axes of \ to the orientation of the particle is same for every particle. 
The total magnetic field H(r) in (E7J) is 

H(r) = IP(r) +H*(r), (48) 

where H p (r) is the field from permanent magnetization density M p , and is given by ([!]) with M set equal to M p , 
whereas H'(r) is due to the induced magnetization: in (|l|) set M equal to M 1 . Note that even an isolated particle 
will have an induced magnetization because the demagnetizing effect will give rise to a non-zero H p , and the total H 



inside the particle must be determined self-consistently, as it both induces a magnetization, (47), and is partly (H l ) 
determined by that magnetization. 

Because of this requirement of consistency between total field and magnetization, the magnetic interaction of 
polarizable particles is a multi-body interaction, and is much more complicated than the pairwise multipole interactions 



of permanently magnetized particles discussed in section IV A. To write down the interaction energy Ti M of a 
configuration of polarizable particles, it is convenient to calculate the work done assembling the configuration, starting 
with the particles well separated from each other at infinity. For any arrangement of particles, the total magnetic 
energy is 

H T = — [ d 3 r |H p | 2 - - [ d 3 r W • M\ (49) 
8tt J 2 J 

The first term is the demagnetization energy (|^) of the permanent magnetization distribution, and the second term 
represents the work done in introducing linearly polarizable material into this permanent field. Evaluating (^9|) for 
an isolated particle defines the self energy E sel f per particle. The difference 

H M = H T - NE self (50) 

equals the work done bringing initially isolated polarizable particles into interaction with each other. 

The stability @ and hence the lower bound @ for this system follows from the positivity of field energy. Rewrite 
the magnetic energy in equation (B9h as (see appendix A for details) 
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H T = J d 3 r H • (1 + A-kx) ■ H. 
Since 1 + Airx is positive definite, H T > 0, hence 



H M >-NE 



self 



(51) 



(52) 



The stability and lower bound then follow as in the case of identical hard core particles. 

To p rove the upper bound (|l0|), consider two subsystems 1 and 2 separated by a distance d > d > as in 
section 



III A[ Write the interaction Hamiltonian between these two subsystems as 



n/M 



H{ 



M 



(53) 



where H\ l and Hif are the magnetic interaction Hamiltonians for the subsystems by themselves; i.e. Hf 1 is the 
energy of subsystem 1 with subsystem 2 placed infinitely far away. We show in appendix B that, for positive-definite 
X, 



n 12 = u 1 



M 



■ATX, 



(54) 



where is odd under reversal of the permanent magnetization of particles in subsystem 1 by the 9 operator, and 
M12 < 0. Since Off is odd, its ensemble average vanishes. Then 



<«&>a=o = <A/&)a=o<0 



(55) 



which is sufficient to prove the upper bound on F, as discussed in section [II C. 



D. Quantum Systems 

Our proofs extend to quantum systems. Consider a system of N identical spin S particles in volume V, which may 
obey Boltzmann, Fcrmi-Dirac or Bose-Einstein statistics. The Hamiltonian is 



H = K + H 



n 



EX 



H 1 



(56) 



Here K = Y^iLiPi/^m is the kinetic energy operator. The exchange interaction ||l| is H EX = \ J( r ij)^i ' Sj, 
where Si is the spin operator for particle i, rij is the distance between particles i and j, the couplings J(fij) are 
assumed to satisfy conditions consistent with stability, and the sum is over all i ^ j. The dipole interaction H M 
is given by (27) with dipole moments /J, i = gSi. The residual interaction H R may be any stable and tempered 
interaction that remains unchanged under simultaneous spin reversal of all particles, such as the hard core or central 
force interactions discussed in sections IV A and IV B . Implicit in our definition of the Hamiltonian is the confinement 
of particles in a region 1Z of volume V with hard wall boundary conditions. The partition function is 



z = I Tr{f 



-H/k B T 



} 



(57) 



For 



where the trace Tr is carried out over states of appropriate symmetry with respect to interchange of particles. 
Boltzmann particles C = N\ and for fermions and bosons C = 1. The free energy F is —ksTlogZ. 

The stability of H, in the sense that its spectrum has a lower bound proportional to N, see ([13]), follows from the 
positivity of K. (no negative eigenvalues) and the stability of H M + H EX + H R . Lower bounds on classical energies of 
the sort derived in sections IV A and IV B| are easily extended to operator inequalities which demonstrate the stability 
of H M + H R . That this stability is preserved upon adding H EX requires a suitable choice of J(r). For example, in 
the case of hard core particles it suffices that J(r) be bounded and decrease more rapidly than r _3_c for some e > 0. 

We now address the upper bound on the free energy. Write the Hamiltonian in the form 



H — Hi + H2 + H12 



(58) 



where Hi includes all terms in ( |56|) (one particle, two-particle, etc.) involving only particles in the set Si with labels 
1, 2, Ni, H2 includes the terms involving only particles in the set £2 with labels N\ + l,N\ + 2, N\ + N2, and 
H12 all the remaining terms. The interaction energy H12 is a sum of magnetic dipole, exchange, and residual terms 



H12 = H 



M 



n_/EX 1 n_j R 

At 19 + rt 19 - 



(59) 
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It will also be convenient to define a scaled Hamiltonian 

H(\) = Hi+H 2 + \Hu, 0<A<1. (60) 

To begin with, we assume that the Ni particles with labels in Si are confined to a region of space 72i C 72 with 
volume Vi, while the N2 — N — N\ particles with labels in S2 are confined to a region IZ2 C 72. of volume V2, and 
that the minimum separation of 72 1 and 722 is a distance d > do. Note that (in general) 72 is larger than 72 1 U 722. 
The Hilbert space of this system is spanned by product states of the form 

l<m> = I<M1, ^l)) \Xn(Nl + 1, Ni + N 2 )). (61) 

where the integer arguments are particle labels. The {|</> m )} form a complete set of states for the Ni particles in IZi 
with appropriate symmetry under interchange of particles and {|xn)} is a similar set for the N2 particles in IZ2. The 
product states (^) have no symmetry for the interchange of particles between the two sets Si and 1S2 > as indicated 
by the superscript U. Due to this lack of symmetry the partition function is given by 

zZn, = ^Tr { e-^} (62) 

where D — N\\N 2 \ for Boltzmann particles and D — 1 for fermions and bosons. If this partition function is evaluated 
using TL with A = in (|60|), — fc^T times its logarithm is the sum F\ + F2 of free energies for the separate systems of 
particles in regions 72 1 and 1Z 2 , each evaluated as if the other region did not exist, since the interactions between the 
two sets of particles have been "turned off" . 

An upper bound on the free energy F of the full system of N particles confined to region 72, in the form of the 
equation (|l0|) is obtained through the following steps. 

1) "Turn on" the interaction between the particles in regions 72 1 and 722 by letting A increase from to 1. The 
resulting free energy is denoted by F N 7 N . The superscript C indicates that particles arc confined to regions TZi and 

n 2 . 

2) Remove the constraint that only particles in set Si are found in 72 1 and only those in set 1S2 are found in 722. 
Any particle may be anywhere in the union of 72 1 and 72.2, provided there are exactly N\ particles in 72 1 and N2 
particles in 72-2. This involves introducing appropriately symmetrized wavefunctions. 

3) Relax the constraint that precisely N% particles are in Ri and N 2 in 72.2. We still require that there be a total 
of N particles in the union of IZi and 72-2- 

4) Relax the constraint that the particles lie in either 72i or 722 , so that all particles can be anywhere in the larger 
region 72. 

Since steps 2, 3, and 4 do not require special attention to long range interactions, they are discussed in Appendix 
C. Basically, each time a constraint is relaxed the free energy decreases, except for step 2, where it remains constant. 
Thus the upper bound obtained for F^ N applies to F. 

C IA 

For step 1 we use the fact that F N ^ N ^ is a concave function of A, and therefore 

Ff&r a (1) < F C N » N2 (0) + (dF c ^ N2 /d\) A=0 , (63) 
which is known as the Bogoliubov inequality p3| . The first term on the right is F\ + F2 , and the second is 

(Hu)x=a = (Wf 2 f )A=o + {H? 2 X )x=o + (Wf 2 )A=o, (64) 
where the averages are with respect to Boltzmann weights with A = 0. For example 

(H M 2)x=o = Tr e -(Wi+w»)/*BT }/Tr { e -(.-H 1+ n a )/k B T } (65) 

The third term in ( |64| ) is bounded by the upper bound A 12 on TC^- The first and second terms vanish for the following 
reason. Let 9 be the anti-unitary spin reversal operator which reverses all the spins of the particles in collection S\. 
It is a symmetry of Tt\, because Tif 1 and 7if x involve products of two spin operators and Tlf is, by assumption, 
invariant under the reversal of all spins. Consequently, if the {|</>m}} are the eigenstates of 7Yi, the states 9\<j) m ) = \4> m ) 
are also eigenstates with the same eigenvalues. In evaluating the trace in the numerator of ( |65[ ) we can employ a 
complete set {|<^m)|Xn)}, where the {|Xn}} are the eigenstates of H2, or, equivalently, {\4>m)\Xn)}- However, 

(4>m\(Xn\n%\4>m)\Xn) = -(<Pm I <Xn|W^ ro ) \Xn) (66) 

because 
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(0 m |S 4 |0 m ) = -(<p m \Si\<f>m) (67) 

for 1 < i < iVi, and TLi 2 1S a sum of pairwise products of spin operators, one from the collection Si and one from £2- 
Thus (7ij|}A=o is equal to its negative and vanishes. The same argument applies to {Hf 2 x )\ = Q. 

Consequently, after following the remaining steps 2, 3, and 4 contained in appendix C, we have a bound corre- 



sponding to (10) in the classical case. This bound together with that in (H), which (as already noted) follows fro m the 



stability of H, completes the proof of the existence of a thermodynamic limit in the quantum case, see section [II A . 

It is possible to treat the translation degrees of freedom classically and the spin quantum mechanically. The proof 
of a thermodynamic limit for such a "semi quantum" model is similar to the "fully quantum" treatment. The averages 
of Tif 2 and TCf 2 x involve sums over all spin states and integrals over all particle positions. Using the spin reversal 
operator 9, introduced earlier in this section, and doing the sums over all spin states first, one sees that the average 
values of Hi 2 and Tif 2 x vanish. The residual interactions T-L p 2 are bounded by A12, and therefore our proof goes 
through. 

One may define another "semi quantum" model of classical spins and classical dipoles with center-of-mass motion of 
the particles treated quantum mechanically. The averages of Hff and 7if 2 x contain integrals over all spin orientations 
and sums over all spatial wavefunctions in this case. The proof is similar to that of classical particles. Doing the 
integrals over particle spins first, gives zero for the average of 7i\ 2 + Hf x , because they are both odd with respect to 
spin reversal. The sum over the spatial wavefunctions then yields t he des ired upper bound on (TL\ 2 )\ = q and our proof 
goes through. We can even include classical polarizability (section | TV C ) in such a model. In this case the average of 



TLi 2 is non positive, and again we get the required upper bound (|10|) on the free energy. 

V. ELECTRIC POLARIZATION 

So far we phrased our discussion entirely in terms of magnetic interactions. However our proof applies to many 
electrically polarized or polarizable materials as well. Electric fields E and D fields satisfy the same Maxwell equations 
as magnetic fields H and B, respectively, provided no free charges and currents are present. By replacing fields H, 
B and magnetization M with fields E, D and polarization P, respectively, our proofs run exactly as for magnetic 
materials. We assume stability, tempering of the residual interactions, and the existence of a 9 operator. Thus we 
prove the existence of a shape independent thermodynamic limit for the electric analogue of each of the classical 
models discussed in sections |y| 

For ferroelectric identical hard core particles, use the Hamiltonian 

U = H P + H HC (68) 



where TL P is the electric interaction between particles. The proof follows section [IV A with the replacement of fields 



discussed at the beginning of this section. For electric dipolar systems with central forces the interaction Hamiltonian 

H = H P + H {n) + H central . (69) 



The proof goes through as in section |IV B[ with the existence of a shape independent thermodynamic limit for electric 
dipolar soft spheres and Stockmayer fluids as special cases. For electrically polarizable particles, H = Ti p + H HC 
where TL P , the electric interaction between polarizable particles, is defined analogously with H M in equation (|5C|). 
After the repla cements of fields and magnetization discussed at the beg innin g of this section, the proof proceeds as 
in section IV C . Quantum models can be treated as discussed in section [VP . 

Adding another layer of complexity, we may define models combining features of electric and magnetic models 
already discussed. The possible variations are too numerous to describe individually. We simply note here that 
stability including only magnetic or only electric interactions ensures stability with the two added together. Finding 
a 9 operator may be more difficult. For example, a 180° rotation axis must lie perpendicular to both M and P. 

For the crude dipolar particle models discussed above, the electric and magnetic problems are fully equivalent. For 
applications to realistic models of specific materials, however, it is generally harder to find a 9 operator in the case 
of electric materials. The analogue of time reversal (which provides a 9 operator for superparamagnetic particles) 
is charge conjugation. This has the undesirable effect of turning matter into antimatter! Without a 9 operator we 
cannot apply our proof. 

For an example of a model outside our proof, consider a fluid of H 2 molecules (water). Modeled as a dipolar 
hard sphere or a Stockmayer fluid the thermodynamic limit follows from our discussions above using rotation as the 9 
operator. Real H2O molecules lack this rotational symmetry. For example, the TIPS 3 site model of water p3] places 
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positive charges on the hydrogens and a negative charge on oxygen. Coulomb interactions between all intermolecular 
pairs of charges, and a Lennard- Jones interaction between the oxygen atoms, gives the interaction 



Ti-mn = + ^2 6 ' (^) 

ij ■> 



between any two water molecules m,n, where i and j run over the charges qi and qj, respectively, on molecule m 
and n, rij is the separation between charges and roo is the separation between oxygen ions. For the above model 
there is no symmetry operation available which can reverse all the electric interactions while keeping the residual 
(Lennard-Jones) interactions unchanged. A 9 operator therefore does not exist for this model and we cannot apply 



our proof. We suppose that our qualitative argument (section [I B ) demonstrating existence of a shape independent 
thermodynamic limit based on domain formation still holds, but for technical reasons we cannot prove it. 

Yet another example is provided by models of ferroelectric materials with mobile charges. Charges transfer among 
dissimilar chemical elements, so no 9 operator is likely to exist. Furthermore, the definition of polarization becomes 
ambiguous, depending on how the unit cell is defined in the case of a crystal |2q ], or on the surface charge for non- 
crystalline materials. The question of appropriate thermodynamic limits for model ferroelectrics with mobile charges 
remains under discussion fl26| . 

For models containing only bound charges, our qualitative argument (section [II Bj ) suggests the equilibrium state 
has no depolarizing energy density. When the microscopic coulomb charges are taken into account, including the 
possibility of molecular dissociation and the Fermi statistics of the constituent particles, the problem falls into the 
class of materials for which Lieb |^7j proved a thermodynamic limit. Free electric charges screen the 1/r Coulomb 
potential, restoring the thermodynamic limit. 



VI. CONCLUSIONS 



We proved the existence and shape independence of the free energy density for a variety of dipolar systems. Three 
essential conditions were identified: stability, tempering of residual interactions, and the existence of a 9 operator 
that commutes with Tit, TL2 and 7if 2 while reversing the sign of Wff . Our proof covers systems of identical hard core 
particles with uniform permanent magnetization. We also treat dipolar soft spheres and Stockmayer fluids, systems 
with magnetizable or polarizable material, and we consider electric as well as magnetic dipoles. Except for the case 
of super-paramagnetic particles, the existence of a 9 operator requires symmetries such as a 2-fold axis of rotational 
symmetry perpendicular to the magnetization/polarization of each particle. 

Having proven shape independence of the free energy we now consider some implications of the proof. For the 
systems covered by our proof, thermodynamic states and phase diagrams do not depend on size or shape [ p8[ . 
Intrinsic thermodynamic quantities such as pressure and chemical potential are independent of sample shape and 
position within a sample. When calculating free energies of magnetized states, care must be taken to remove the 
depolarizing field if uniform magnetization is assumed. Failure to do so leads to either boundary-condition or shape 
dependence of thermodynamic properties p^j2^ -|30[| . Two convenient ways to remove the demagnetizing fields are to 
study highly prolate ellipsoids or to use tin- foil boundary conditions ||. 

Computer simulations []3l| suggest dipolar fluids such as ferrofluids may spontaneously magnetize. Experiments 
on supercooled CoPd alloys ^] claim to observe a spontaneously magnetized metastable state. How then can we 
reconcile shape independence of the free energy with these reports of spontaneously polarized liquid states? Domains 
or textures must form, as in equilibrium solid ferromagnets. Liquids lack the crystalline anisotropy required to sustain 
a sharp domain wall (see Fig. lb). De Gennes and Pincus 33 note, however, that domain wall thickness should 
be comparable to system dimensions in a magnetic liquid. We conclude that a spontaneously polarized liquid has 
a position dependent axis of polarization that rotates so as to lie tangent to the sample surface. Since vector fields 
tangent to the surface of simply connected volumes exhibit singularities (the "hairy billiard ball" problem), magnetized 
fluid droplets must contain defects in the magnetization field M. Possible textures include line (Fig. lc) or point 
defects (not shown). Away from such defects the magnitude of magnetization |M| is independent of position within 
the sample. In computer simulations a uniformly polarized state is observed because the Ewald summations || drop 
the surface pole energy (via "tin-foil" boundaries) and thus mimic an infinite, boundary-free, medium. 

Our proof is valid only in zero applied field. The case of an applied field is still open. The free energy in a field 
depends on shape as outlined in section [n]. The free energy increase due to the demagnetizing field causes a droplet 
of paramagnetic liquid to elongate in an external field, minimizing its magnetic energy by reducing D [pij. However, 
we conjecture the existence of an intrinsic thermodynamic energy, defined in equation (J5|) by subtracting the shape 
dependent demagnetizing energy from the full shape dependent free energy. 
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Finally, let's consider some dipolar systems to which we are unable to apply our proof due to the lack of a 6 
operator. Particles with non-symmetric shapes and magnetization fixed with respect their body (Fig. |^a) provide an 
example. Each particle is a cube with protrusions (conical, hemispherical and cubic) on three faces and matching 
indentations on the opposite three faces. Each face of one particle fits exactly into the corresponding opposite face 
of another particle. Recall that we exploit symmetries of the Hamiltonian when applying a 9 operator to the internal 
coordinates 0$. The 6 operator, if it exists, reverses the sign of Hf^ while leaving 7ii invariant. Rotations are not 
a symmetry of the internal coordinates for these particles. There is no evident symmetry of the Hamiltonian which 
could be used as a 8 operator. 

When these particles are tightly packed (i.e. the limit of infinite pressure and hard cores), they align parallel to each 
other. M is uniform throughout space and therefore the thermodynamic limit does not exist. At finite pressure it is 
possible to accumulate enough free volume to insert domain walls, (Fig. |^b), so we expect that the thermodynamic 
limit does exist even though the conditions of our proof are not obeyed. That is, our conditions are sufficient but not 
in all cases necessary. 



Another interesting example is provided by hydrogen absorption in metals 35 1 . Interstitial hydrogen creates elastic 
strain fields that fall off as l/r 3 , and hence are referred to as dipoles. This is a rather confusing notation, since 
the angular dependence is actually a sum of a quadrupole and a monopole contribution. The monopole term is the 
trace of the strain tensor and governs lattice dilation. Provided the system remains coherent (the lattice structure 
is dislocation free), the monopole interactions are shape dependent, attractive and infinite ranged. The interactions 
reduce the energy but create no force. In the coherent state, the phase diagram for liquid-gas transitions of hydrogen 
in a metal depends on the shape of the metal Q . 

The coherent state itself is a metastable state. The true equilibrium state is incoherent, with dislocations relaxing 
the lattice strain. If one allows for dislocations, the elastic interactions become short ranged and the thermodynamic 
limit is restored. The sample itself, however, may have disintegrated into a fine powder! 
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APPENDIX A: TOTAL MAGNETIC ENERGY FOR POL ARIZ ABLE DIPOLES 

We start with equation ( pi] ) for TL T and show that it equals equation ([l9]). Use equations p?| ) for induced magne- 
tization M l and (^) defining the permanent and induced fields H p and EP to rewrite equation ( p3l| ) as 

H T = — [ d 3 r |H p | 2 + — [ d 3 r IPTI 2 + — [ d 3 r H p ■ H 1 (Al) 
+ 1 [ d 3 r W ■ M l + I f d 3 r U l • 1VP 



For any two magnetization distributions M^r), and M 2 (r) and the fields Hi(r) and H 2 (r) caused by them respec- 
tively, the following identity holds B : 



i- J d 3 v Hi H 2 = 1 J d 3 v H!-M 2 = -i J rf 3 rH r M! 



(A2) 



Since equation ( |A2| ) holds for any t wo arbitrary magnetization distributions, we set Mi = M 2 equal to the induced 
magnetization M 1 . Then equation ( [A2] ) gives 

— / d 3 r |H ,: | 2 =-- I d 3 r ■ M\ (A3) 



8tt J 1 2 _ 

Similarly, setting Mi = M p and M 2 = M 4 in equation (A2) gives 



-!- [ d 3 r H p • H 4 = -I I d 3 r H p ■ M\ (A4) 

87T ./ 2 



Usin g eq uations (A3) and (A4) to simplify equation (Al) gives equation (|49"|), proving equality of our expressions (|4£ 
and (pi]) for TL T . 
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APPENDIX B: INTERACTION ENERGY BETWEEN TWO SUBSYSTEMS OF POL ARIZ ABLE 

DIPOLES 



Let Hj(r) and H^r) be the fields due to the permanent polarization in subsystems 1 and 2 located in non- 
overlapping regions IZi and 7^-2 • The induced magnetizations in the two subsystems can be written in the form 



Mi(r)=M;(r)+M 1 (r), 
M^(r)=M^(r)+M;(r), 



(Bl) 



where Mf(r) is the induced magnetization which would be present in subsystem 1 were subsystem 2 absent, and 
M^r) that of subsystem 2 were subsystem 1 absent. 

for TL T and ( pp[ ) for 7i M to find the interaction Hamiltonians Hf 1 , Ti.% 1 and Ti. M for the two 
role system, respectively, we write the interaction energy in (p3h as 



Using equations 
subsystems and the w 



njM 

H-12 - 



J d 3 r H? • Hf - i J d 3 r H p 2 ■ M* - i J d 3 r • M s 2 

-\ J dh (h? + h?) • m; - \ J d 3 v (h? + Hf ) • m 2 . 



(B2) 



Break Tif^ into odd and non-positive components Hf^ — Ofi + i where 



AT 



4-7T 



d 3 r • M? 



d 3 r H? • Mo 



is odd under reversal of the permanent magnetization of particles in subsystem 1 by the operator, and 



- ,v 12 — 



\ Jd 3 r (H?+H^)-M;-i|d 3 r (Hf + Hf ) • M 2 



(B3) 



(B4) 



is non-positive as we now show. A theorem by Brown Q states that for a paramagnetic polarizable material in 
an applied field, the unique induced magnetization M' given by equation ( f47| ) minimizes the total magnetic energy. 
Applying that theorem to our system we observe that 



\ J d 3 r (H? + H|) • Mi - I J d 3 v (H? + Hf) • M 2 < 
-i | d 3 r (H? + H*) • Ml - \ J d 3 r (H? + H?) • MJ. 



(B5) 



because the interaction-induced magnetization M 1 has lower energy than the isolated self-magnetization M s . Upon 
replacing M^ in (p4|), with Mj - Mj, (pl|), one sees that ( p35| ) implies that 



< 0. 



(B6) 



APPENDIX C: QUANTUM SYSTEMS 



Here are the details for steps 2, 3, and 4 in section IVD. For step 2 note that the formal Hamiltonian 7i(A = 1), 



defined in equation (|60|), is symmetrical under the interchange of any two particles, since the right side of (|5§|) is 
simply a way of segregating terms in the sum representing TL. To allow any particle to be anywhere in 1Z\ U 7^.2, 
subject to the requirement of N\ particles in 1Z\ and N2 in IZ2, replace the Hilbert space of the form (^) with another 
spanned by states with appropriate symmetries under interchanging of any pair of particles. We now construct such 
a Hilbert space for each type of statistics as indicated by Fisher || . 

First consider identical particles obeying Boltzmann statistics and let j=l,2... be a complete orthonormal 

set of single particle states (including spin) for a particle confined to 1Z\. A basis {|</>m)} for the particles with labels 
in 1S1 can then be written in the form 

1^(1,...,^)) = | A1 ™ 1 (l))l^ m2 (2))lM m3 (3))... (CI) 
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where m stands for the sequence (mi,m 2 ...m]y 1 ) of integer labels. In the same way, a basis {|Xn)} f° r particles with 
labels in S 2 can be constructed using single particle states {|^ fe )}, k=l,2... for a particle confined to 1Z 2 . 

To construct a basis in which any Ni particles are in IZi and any N 2 particles in 1Z 2 , we proceed as follows. Consider 
the collection V of permutations p of the integers (1, 2...N), where p(j) is the image of j under p, with the property 
that 

p(l) < p(2) < ... < p(Nx), (C2) 
p(Ni + 1) < p(Ni + 2) < ... < p(Ni + N 2 ). 

It is clear that V contains P — N\/(Ni\N 2 \) permutations, one for each way of partitioning the integers from 1 to N 
into two collections, one containing N\ and the other containing N 2 integers. Then define states 

\lpm,n,p) = CT p \<f> m {l,...,N 1 )) \x n (Ni + l,...,Ni+N 2 )), (C3) 

where the operator a p applies permutation p to the N\ + N 2 arguments. The set {\ip m ,n,p)} for m and n defined 
previously, and p belonging to V, forms an orthonormal basis for the N particles in IZi U 1Z 2 , allowing any particle to 
be in either region, subject to the constraint of N± particles in TZi and the remaining N 2 in 7Z 2 . 
To prove orthonormality, 



note that factors 6 mm >,5 nn i follow from orthonormality of the single particle states in (CI). To get the factor S pp i 
consider p p' bot h belonging to V '. There is at least one particle I which is in TZi under p, and in TZ 2 under p' . The 
inner product (C4) contains a factor (^ mi (p(l))\v ni (p'(l))) that vanishes because the vanish outside TZi, the \v k ) 
vanish outside 1Z 2 . Recall that TZi and 7Z 2 do not overlap. Finally, the normalization condition in (C4) follows from 
the unitarity of a p . In addition note that for p ^ p' both belonging to V, 

(lp m ,n,p\'H% u N 2 \'4>-m>,n>,p>) = (Xn I (0m Wp 1 H% 1 t N 2 a p' l^m' > \Xn' ) = 0, (C5) 

because the Hamiltonian does not interchange particles between the two regions TZi and 7Z 2 . 
We can use the set of states {\ipm,n, P }} to evaluate the trace in the partition function (|5~ 



Z k,N 2 = ^EE^-I^K" 1 e~ H ^ k * T a p \M\ Xn ). (C6) 

m,n ptP 

The Hamiltonian, being symmetric, commutes with a p . The sum over p therefore just gives a factor of P, so that 

Z° NuN2 = j^jEfrnlfoml e- n ^' kBT \4> m )\Xn), (C7) 
^ ^ m,n 

which is the same as the partition function Z^ N ^ (defined in equation (|62|)) evaluated in the Hilbert space spanned 

by {\^m,n)}- The free ener gy F N t ,N 2 therefore is equal to F^ N2 - 

For fermions (— ) and bosons (+) the appropriately symmetrized orthonormal states are 

K„}=^^E( ± ) 7r(rt l^»( 1 '-' Ar i)) \Xn(Ni + l,...,Ni + N 2 )), (C8) 

peV 

where ir(p) is for an even and 1 for an odd permutation p, and \4> m ) and \\ n ) are assumed to have appropriate 
symmetry with respect to interchange of particles within TZi and within 1Z 2 , respectively. We use the set of states 
{IV'm «)} to evaluate the trace in the partition function (|57j): 

Z Ln 2 = J2jsY,J2 (±r (P) W) (Xn I <0 m I V 1 Vp\4>m)\Xn). (C9) 

m,n peV peV 



Only terms with p = p' survive because of (C5). Since (±) 2?t (p) = 1, the rest of the proof is similar to the Boltzmann 
case. 

Step 3 in section IV D follows from the observation that the partition function Z for N particles in TZi U 1Z 2 is a 
sum of positive terms of the form 
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N 

z ° = E z k,N- Nl , (cio) 

Nl=0 

where Z%j N ^ is the partition function for Ni particles in IZi and N2 particles in IZ2 ■ Consequently Z c is not smaller 
than Z < f ii N2 for any particular Aq and N2 = N — Ni , and 

F c <F% uNr (Cll) 

Step 4 in section IV D| exploits the fact that region 1Z contains IZi U TI2 ■ Whenever particles are allowed to move in 
a larger region, the partition function goes up and the free energy goes down. The quantum version of this result is 
based on the fact that an energy eigenstate \ip m ) of 7i in the smaller region, with energy E m , belongs to the Hilbert 
space of the larger region. Although it is not an eigenstate of 7i in the larger region, it is still the case that 

^mWm) = E m . (C12) 

Now we make use of Peierl's theorem, which states that 

Tr e -H/k B T > e -{<t>lYH/k B T\4>i) ( C13 ) 

where {|0z}} is any orthonormal set of states. Choose to be the set of energy eigenstates {|VVn)} of 7i in 

1Zi U 7^2 ■ It follows that the partition function Z for region 1Z, is greater than Z c for lZi U7l2, and thus 

F < F c . (C14) 
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FIG. 2. (a) Permanently magnetized particle lacking a 6 operator. The arrows indicate the direction of magnetization, (b) 
Domain formation requires gaps between particles. A cross-section is shown. 
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